More about the Tolman-Oppenheimer-Volkoff equations for the generalized Chaplygin 

gas 
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We investigate the Tolman-Oppenheimer-Volkoff equations for the generalized Chaplygin gas with 
the aim of extending the findings of V. Gorini, U. Moschella, A. Y. Kamenshchik, V. Pasquier, and 
A. A. Starobinsky [Phys. Rev. D 78, 064064 (2008)]. We study both the standard case, where we 
reproduce some previous results, and the phantom case. In the phantom case we show that even 
a superluminal group velocity arising for a > 1 cannot prevent the divergence of the pressure at a 
finite radial distance. Finally, we investigate how a modification of the generalized Chaplygin gas 
equation of state, required by causality arguments at densities very close to A, affects the results 
found so far. 
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I. INTRODUCTION 

The presence of dark energy (DE) component in our 
universe seems to be a matter of fact. In this situa- 
tion, the study of spherically symmetric solutions of the 
Einstein equations with additional matter or geometric 
terms describing DE turns out to be an important issue. 
Moreover, since the possibility of a phantom DE, namely 
a dark component which has an equation of state param- 
eter w < —1, is not completely ruled out by observation 
BS&SH, this study is even more impelling from the 
standpoint of wormholelike solutions. 

In this paper we generalize the results found in [l| to 
the case of a DE described by the generalized Chaply- 
gin gas (gCg) @ (see also [§j]), in particular, in the light 
of our recent paper 0, where the possibility of a super- 
luminal regime has been noticed; in 0] we have shown 
that the superluminal regime does not contradict causal- 
ity provided some modifications of the gCg equation of 
state are made which do not affect the present and past 
evolution of the universe. It is also interesting to investi- 
gate the effects of a very large gCg group velocity on the 
geometries already found in 

In Sec. HH we write down the Tolman-Oppenheimer- 
Volkoff (TOV) equations, and in Sec. IHII we exhibit their 
special constant pressure solutions. Sections IIVI and [V] 
are devoted to the analysis of the normal case \p\ < p 
and of the phantom one \p\ > p, respectively. In Sec. IVII 
we address the superluminality issue and we investigate 
how a modification of the gCg equation of state changes 
the pressure solution without spoiling the results found. 
Section IVIII is devoted to discussion and conclusions; in 
particular, we emphasize the differences arising in the 
gCg model compared to the standard Cg (a = 1) case. 



II. TOLMAN-OPPENHEIMER-VOLKOFF 
EQUATIONS IN THE PRESENCE OF THE 
GENERALIZED CHAPLYGIN GAS 

We assume a static spherically symmetric space-time 
geometry (c = 8nG = l): 1 

ds 2 = e v{r) di 2 - e^ r) dr 2 - r 2 (d9 2 + sin 2 6d(f> 2 ) (1) 

and a perfect fluid stress-energy tensor 

T^v = (p + p)u^u v - g^p . (2) 

Einstein's equations, conservation of the energy and the 
boundary condition 

e^ (0) = 1 (3) 

together imply the TOV [if], [H| differential equation 

dp _ (p + p)(M + 4irr 3 p) 
dr 

where 



r(Sirr - 2M) 
p{u) 47TU 2 du 



(4) 



M(r) = 

satisfies the equation 

^ = 4^r 2 p , M(0) = . 
dr 

The equation of state of the gCg 0, HJ is given by 

A a+1 



P= ~- 



(5) 



(6) 



(7) 



1 In the preceding paper [l| we have used a different normalization: 
G = 1. 
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where A and a are positive constants. By plugging ([7]) in 
the TOV equation we get the following first-order system 
of differential equations: 



d\p\ (A^ 1 - \pf +1 ) (M - Airr 3 \ 
dr 
dM 



- Airr 
dr \p\P 



\p\Pr(8irr - 2M) 



(8) 
(9) 



written in terms of the modulus of the pressure and of 
the parameter j3 = I /a. 



III. CONSTANT PRESSURE SOLUTIONS 

One solution with constant pressure is the de Sitter 
universe: 



P 



-p = -A , M(r) 



7 7rAr 



where the de Sitter radius is given by 



(10) 



(11) 



and the curvature does not depend on (3. 

The Einstein static universe (ESU) is another solution: 

p = -A cff EE -3~^A , 

p = A off + p = 3»+! A , p = 2A oS , (i2) 

M (r) = AttA cS r 3 . 
The radius of the spatial spherical section of the ESU is 




(13) 



In the limit (3 — > oo the two radii te and rds coincide 
(but not the geometries), while in the ultrasuper luminal 
(3 -> limit r E -> y/T/A . 



IV. THE NON-PHANTOM CASE p > \p\ 

Consider the system (JS])-© and assume the dominant 
energy condition p > \p\. 

We fix the pressure p{rb) and the mass M(ri,) at some 
value r = rj,. The following inequality holds: 



Ann > M(r b ) 



(14) 



At these initial conditions the pressure can neither vanish 
nor attain the value p = —A. The former case is excluded 
because, whenp — > 0, the right-hand side (rhs) of Eq. 
is positive and therefore the pressure cannot vanish. The 
value p = —A is also excluded, which can be seen by 
rewriting Eq. ([8|) as follows: 

din (br* _ A*») = -Q + D ^^ dr ■ (15) 



If the pressure could approach —A, the left-hand side 
(lhs) would diverge logarithmically while the rhs would 
stay regular: a contradiction. 

If the condition p > \p\ is satisfied, then p(r) > A or 
equivalently 



A < p(r) < 



(16) 



Equation @ then implies that M{r) grows at least as 
fast as r 3 and, therefore, the equality Airr = M(r ) is 
achieved at a certain radius r$. 

At r — ro the denominator on the rhs of Eq. (JS]) di- 
verges unless rgp(r ) = — 1. 

To elaborate on this point let us rewrite Eq. fSJ by 
expanding the relevant quantities in the neighborhood of 
ro. Let pq and po be the values of the pressure and of the 
energy density at vq. To first order the mass is given by 



M{r) 



Anrn — Anr, 



(17) 



where we have set e = tq 
have 



IpoI' ' 
r. As for the pressure, we 



\p(r)\ = \po\+P(e) , 



(18) 



where p(e) vanishes when e — > 0. 

Hence, for e ~ Eq. © takes the following approxi- 
mate form: 



g = (AW-M^) (1-rgbol) 
de 2e(|p |/3-r2A/3+i) 

A logarithmic divergence is manifest unless 



\po\ 



(19) 



(20) 



The inequality —A < p{r) < provides a lower bound 
for Tq: 



A <p Q 



1 



r > 



(21) 



the inequality p(r) > A implies a bound for the mass 
function: 



M(r) > -7rr 3 A 
Since M(r ) = Atttq it follows that 



<r < 



(22) 



(23) 



[see Eq. (|13[) ]. The approximate equation for p is written 
where 



2rt 



(25) 
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so that 



p = A^fl + 2C e 



(26) 



where A is an integration constant. The solutions 
p(r), M(r) are characterized by the parameters A and tq 
which are in turn determined by the boundary conditions 
at r — Tb- 

In order to avoid the apparent singularity on the rhs 
of Eq. ([5]) it is useful to change the radial coordinate as 
follows: 

r = r sinx- (27) 
The gCg TOV equations take the following form: 
d\p\ _ (A0+ 1 - \pf +1 ) (M - 4 7 rr 3 ,|p| sin 3 x) cos X 



dx 



\p\P sin x(87rro sin x ~ 2M) 



dM , 2 A^ +1 

— = 47rr sin xcosx^ 



(28) 
(29) 



We investigate these equations in the vicinity of the equa- 
tor (x= f ) by the substitution 



(30) 



where 6 is positive. Once more, we expand the pressure 
and the mass, 



\ P \ = \+p(5) 



M = 47rro - 27rr (rgA)' 3+1 S 2 



(31) 



(32) 



to lowest order in S. We can then write an approximate 
equation for p(5): 



d P _ p , r x 



where 



2ri 



(r^r 1 - 3 



The solution of 



p = B5 + dS 2 



(33) 



(34) 



(35) 



where -B is an integration constant. The presence of 
the parameter (3 does not appreciably affect the equa- 
tor crossing being of subleading order in 8. The Ein- 
stein static solution corresponds to the choice B = and 
d =0: 



37J+T 

~1T 



(36) 



The solution (f35|) can now be continued to negative values 
of 6 and can describe the equator crossing. This can be 
done by introducing the variable 



1 



sin x 



so that y 6 [l,oo). Equations (|2"51) - (f2"9")) become 
d\p\ (A' 3 + 1 -b|^+ 1 )(My 3 -47rr 3 b|) 



dy 
dM 
dy 



|p| V(8tito - 2My) 



4^r|A^+i 



(37) 

> (38) 
(39) 



Following [l[ we now show that the quantity 87rro — 2My 
is always positive. 

The rhs of Eq. (f39|) attains the smallest possible value 
at p = —A; this implies that 



87rr -2My^ 87rr - 



iirAr 



The rhs is positive for 



y 2 > 



(40) 

(41) 
and (1411 



Since by definition y 2 > 1, the inequalities 
imply that 87rro — 2My is always positive. 

We are now in a position to investigate solutions of the 
TOV equations for X ~ * 7r ; i- e - V ~ * °°- 

Case 1. Consider solutions for which the pressure at- 
tains, at infinity, a generic value in the range < \poo\ < 
A. For such a solution to exist, the integral 



y , Mv 3 - 47^3, |p I 
dv- 



///, 



v 3 (8irr - 2Mv) ' 



(42) 



must converge for y — > oo, where yb = ?'o/ r fa; cannot 
be either positive (otherwise 87rro — 2Mv would become 
negative) or negative (otherwise the integral would di- 
verge logarithmically). Therefore, must necessarily 
vanish. Then, Eq. (|3"9"|) gives the asymptotic behavior 



M 



4^r|A^+i 
3y 3 boo|' 3 



(43) 



and, from Eq. (J3S 



_ 2 (A^-| Poo r i )(A^+ i ~3b 0O r i ) 

dy ~ T ° 6|i>oo| a V 

(44) 

Case 2. \p{y a )\ = at a certain y a > 1. Let M(y a ) = 
M a . In the neighborhood of y a Eq. (|3"8)) can be written 
as 



dy 



A /3+1 M„ 



|p|0(8tito - 2M a y a ) 



(45) 



In order to have a nonpositive derivative it is necessary 
that M a ^ 0. First consider the case M a ^ 0. Assume 
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the following asymptotic behavior for the pressure in the 
neighborhood of y a : 



\p\ ~ D(y a - y)i , 

where D and 7 are positive constants. Equation 
gives 



(46) 



DliVa ~ V) 



and therefore 



7-1 



Df 3 (y a -y)^(87rr -2M a y a ) 



1 



' 0+1 ' 

D/3+1 = ( / 3 + l)A^+ 1 M„ 
87rr - 2M a2/a 

If M a = we assume near y a : 

M(y)~M 1 (y a -yf , 

\p\ ~ £(y Q - y)" , 



(47) 



(48) 



(49) 



(50) 



(51) 



where Mi and -E are positive constants and < £ < a. 2 
Substituting l[50" |) -(|5T ]) into Eq. ^ we get 



Mi((y a - y) 



c _ x _ 4 7 rrgA' 3 + 1 



yiEP{y a - y y° 



It follows that 



and 



„ 4 7 rr2A' 3 + 1 
Mi^p = °- A 



C = l-0a. 
Since £ > we can infer that 

1 1 

JTT <a< P 



Plugging (J50J) and (JSTjl into Eq. (J45J) one finds 

A /3+1 Mi 



Therefore, it follows that 



i + C 

13 + 1 



(52) 



(53) 



(54) 



(55) 



(56) 



(57) 



2 This inequality comes from Eqs. J38P - (130 I t written as 

d\p\ _ (AP+ 1 - \p[ 3+1 ) (My 3 - 47rrg|p|) y 
dM ~ 47i-r-3A/5+ 1 (87rr - 2My) 

which vanishes for y — » y a . 



and 

Ef3+1 = AWM, ^ 
87rrocr 

Combining (|57 )) -([5g ]) with |[53" ]) -([54" P we obtain 

2 



(58) 



and 



1 + 2/3 
1 

1 + 2/3 
167rr 



(59) 



1 A' 3 + 1 (l + 2/3) ' 
2/3+1 _ ^AW)(l + 2/ 3) 2 
4y4 



(60) 



Case 3. p^ = —A. In this case, rewrite the pressure 
as follows: 

H=A-|p|. (61) 
Equation (|38p has the following asymptotic form: 
d|ji| C9+l)|p| 



2/y 



whose solution is 



y 2 



(62) 



(63) 



where F is a positive integration constant. 

V. THE PHANTOM CASE p < \p\ 

The violation of the dominant energy condition is in- 
teresting for the following two reasons: 

1. p < \p\ is a necessary condition for the existence of 
wormholes solutio ns Jl3 . [l3| (but it is not sufficient 
for their stability [l4j). 

2. The possibility of a phantom DE has not been ruled 
out by observation 0, d H i, @| . 

We consider Eqs. ©-© together with the assumption 
p(r b ) < -A . (64) 

Case A. If 



M(r b ) - 4irrl\p(r b )\ < 



(65) 



then d\p\/dr > at r — r b - Since \p\ grows and satisfies 
condition (j64|) . the term M — 47rr 3 |p| stays negative also 
for r > r&. Therefore, d|p|/dr > for r > 7v Then, 
there are three possible asymptotic behaviors: 
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1. For r — » oo, |p| tends to a certain finite value 
|p(oo)| > A. 

2. For r — > oo, |p| diverges. 

3. For r tending to a finite value, say r\, \p\ diverges. 

The first subcase is ruled out since the rhs of Eq. ([5]) 
diverges for r — > oo. 

The second subcase cannot take place because Eq. flSJ) 
has the following asymptotic form: 



dr 



1 



whose solution 



\p\ 



Mp(n)\ 



4+bK)|(r 2 -r 2 ) 



(66) 



(67) 



diverges at the finite distance r 2 = r 2 + 4/|p(rj)|. This 
is a contradiction. 

Then, we are left with subcase 3. In the neighborhood 
of r = n, assume for the pressure the following power 
law behavior: 



\p\ = 



Pi 



(68) 



where r] > 0. Let M(n) = Mi; Eq. flSJ has the following 
form in the neighborhood of r = r*i: 



(n - r 

from which 
and 

Case B. If 



VPi _ +Kr\p\ 

~ 87rri - 2M a 



in 



r)- 2r > 



Pi 



8?rri - 2Mi 



47rr 2 



M(r fe ) - 4?rrf \p(r b )\ > 



(69) 



(70) 



(71) 



(72) 



then d\p\/dr < at r = rj. From Eq. ([5]) it is easily seen 
that |p| cannot decrease ad libitum because at a certain 
finite radius, say r 2 , the term M — 4irr 3 \p\ changes sign, 
becoming negative. 

Then, the pressure has the following profile: fromp(r{,) 
it grows up to a certain (negative) value, say p max at 
r = r%, and then decreases according to the behavior 
described in subcase 3 of A, diverging at r = n with 
n > r 2 . 

In general, p max ^ —A. We prove that the equality is 
ruled out. We rewrite the pressure equation as follows: 



din (\pf +1 - A^ +1 ) = -{p + 1) ^Z^yi dr . ( 73) 



If M - 47rr 3 |p| > and 8nr - 2M > 0, we have a loga- 
rithmic divergence on the lhs while the rhs is regular, i.e. 
a contradiction. 

If we demand that 8irr — 2M vanishes in Eq. (|T3"|) 
has the following form in the neighborhood of r = r 2 : 



dln(\pf +1 -A f3+1 ) 
whose solution is 

| p |/3+l _ A /3+l _ | r _ r2 | 



+ 1 dr 
2 r — r 2 



(74) 



(75) 



r(8nr - 2M) 



For r r2 the lhs vanishes while the rhs diverges, i.e. 
another contradiction. 

The last possibility is that M — 4irr 3 \p\ = and 87rr — 
2M — at r = r 2 - In this case 

P = ~ (76) 



and from Eq. © 

8-Kr - 2M ~ 87r(r - r 2 ) 1 - (r2A) /?+1 . (77) 

Since \p\ > A and r 2 A < 1 the rhs is negative, while 8irr— 
2M > by assumption, which is again a contradiction. 

In conclusion, only two regimes are possible for a stellar 
object immersed in the phantom gCg: 

1. If the boundary conditions at the surface of this 
object satisfy Eq. ([65]) . then \p\ grows and diverges 
at a finite radius r = r±. 

2. If instead inequality (|T2"|) holds true, then the pres- 
sure grows, attaining a maximum (negative) value 
Pmax < — A at r = r 2 , and then it decreases, di- 
verging at r = r\ . 

It is important to stress that, as it can be seen from 
Eqs. ([70| - ([7T|) . the divergence of the pressure at a finite 
radius does not depend on the parameter [3. Therefore, 
the formation of a curvature singularity at a finite value 
of r is unavoidable even for a very large gCg sound speed. 

Then we are in the position to generalize the theorem 
inQ: 

In a static spherically symmetric universe filled with 
the phantom generalized Chaplygin gas, the scalar curva- 
ture becomes singular at some finite value of the radial 
coordinate and the universe is not asymptotically flat. 



A. Wormhole-like solutions for a universe filled 
exclusively with the phantom generalized Chaplygin 

gas 

Let us consider 8-Kri, — 2Mb > and pb < —A and 
investigate the solution of Eq. ||5J) for small values of r. 
There are two possibilities: 

1. 8nr — 2M(r) remains positive up to r = ; 
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2. 8wr — 2M(r) vanishes at a certain r — ro . 

In the first case, let p(r = 0) = po and, using Eqs. |8])- 
@ , we expand the pressure modulus and the mass in the 
neighborhood of r = 0: 



H = \po\ 



12|po| 



2;i 



A/ 



47rA" +1 



3 bo 



/3 



. (78) 
(79) 



In the second case 87rr — 2M vanishes at r = ro. As for 
the nonphantom case of Sec. IIV1 M — 47rr 3 |p| vanishes 
at ro and PoTq = — While in the nonphantom case 
ro is the maximum value of the radial coordinate, in the 
phantom case it is the minimum one and represents the 
radius of a throat. 

We now investigate the crossing of the latter. To this 
purpose, we change coordinates as follows: 



Eqs. 

d\p[ 
dr/ 

dM 
dr] 



r = r cosh 77 . (80) 

become 

(A /3+1 - \pf +x ) (M- 4?rrg |p| cosh 3 7?) sinhr? 



\p\P cosh7y(87rro cosh 77 — 2M) 
A/3+i 



47rr cosh 77 sinh r\ - 



\v\ 



(81) 



The solution of the mass equation near ro, i.e. for small 
rj, is 



M = Airm 



1 



V 



while, rewriting the pressure as 



P 



1 



Eq. (fBTj) takes the following form near the throat: 

d P P , n 
dr] i] 

where 

v/3+l 



Cn 



3~(rgA) ^ 
2rg 



The solution of this equation is 



(82) 



(83) 



(84) 



(85) 



(86) 



where D is an arbitrary constant. It is important to 
emphasize that the wormholelike solutions that we have 
found here are completely different from the Morris- 
Thome- Yurtsever ones fl3( ] since they do not connect two 
asymptotically flat space-time regions, but two regions in 
which the space-time has a singularity. 



VI. THE SUPERLUMINAL NON-PHANTOM 
GENERALIZED CHAPLYGIN GAS AND THE 
CAUSALITY PROBLEM 



The gCg speed of sound is 

, aA a+1 



(87) 



In the range a > 1 (i.e. f3 < 1) it can exceed the speed 
of light for values of the energy density sufficiently close 
to A. This occurs for 



p < a a + 1 A — p s 



(88) 











^max — exp 




+ 1 


w 3.591 











The largest value of p s \ is attained at the maximum of 
the function a°+i, namely at 



(89) 



where W is the Lambert function. At a = a max one has 
/;./ « 1.321A. 

In cosmology a speed of sound exceeding 1 does not 
contradict causality. Indeed, what is required by causal- 
ity is that the signal (i.e. the wavefront) velocity does 
not exceed that of light. As we have shown in Q using 
an underlying tachyonlike field model for the gCg, the 
condition that c s < 1 for the asymptotic "vacuum" state 
p = —p = A is sufficient for the purpose. Therefore, fol- 
lowing [s}, we solve the causality problem for a > 1 by 
smoothly modifying the gCg equation of state to that of 
the standard Chaplygin gas (a — 1) at some value p = p c 
very close to A. 

If we assume that p c < po, where po is the present 
cosmological energy density, the required change in the 
equation of state will occur in our future. However, this 
escamotage does not work in the geometries described in 
the present paper because the change takes place at a 
finite coordinate distance r. 

In order to quantitatively describe the transition in 
the equation of state, assume that p = p c for r = r c and 
suppose that a is now a function of the energy density 
so that a — > 1 sufficiently fast for A < p < p c and it 
becomes constant for p > p c . A simple example of such 
function is 



a 



1 



tanh 



P- Pc 



1 



(90) 



where S is a small (i.e. |A — p c \ 3> S) positive parame- 
ter which tunes the rapidity of the transition and a is a 
constant. 

We now investigate the effect of the transition in the 
equation of state on the pressure and mass configurations 
in order to understand if all the results found so far still 
hold true. 

Assume that p c is close to A, i.e. 



Pc 



A 



A 



< 1 



(91) 
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Then from ||7J) we have 

\p c \=A[l-a(r c )x] , (92) 

to first order in x. Equation (jHJ) now reads 

= (A^)+ 1 -|p|^)+ 1 ) (Af-W| P |) 
dr |p|/ 3 ('')r(87rr - 2M) 1 j 

with a r-dependent parameter /?. We expand all the rel- 
evant quantities in the neighborhood of r c : 

r = r c + e , 

M(r) = M c + 4nr 2 cPc e , 

\p\(r) = \Pc\+P(e), m) 

(3(r) = &+f e. 

Assume that f3 c — 1 and substitute formulas (|9 1 [) - (T92|) 
and the expansions (f9~4"|) in Eq. (|9"5)) . To first order in e 
and a;, we find 

dp 2A (M c - 47rrj?A) (M c - 47rrj?A) 

~de = r c (8^r c - 2M C ) X ~ r c {^r c -2M c ) P ' (95) 

Since d(3/dr does not appear in Eq. ([93)) we can con- 
clude that a rapid variation from the gCg to the stan- 
dard Chaplygin gas equation of state at p — p c does not 
change the results obtained so far. 

We did not consider here the case of the phantom gCg 
(i.e. p < A) because the standard phantom Chaplygin 
gas (a = 1) is already superluminal by itself [see Eq. 
([87]) ]. Moreover, note that, in the phantom case, in order 
that c 2 s < 1 the energy density must satisfy the inequality 

a^TTA<p<A, (96) 

which implies a < 1. 

On the other hand, we have shown in Sec. |V]that the 
pressure has a singularity at a finite radial distance r, 
where the energy density vanishes. Then we conclude 
that in the phantom case a superluminal speed of sound 
cannot be avoided. 



VII. DISCUSSION AND CONCLUSIONS 

In this paper we have investigated solutions of the 
Tolman-Oppenheimer-Volkoff equations for static and 
spherically symmetric configurations of the generalized 
Chaplygin gas, thus generalizing the same analysis car- 
ried out in |l[ for the standard a = 1 case. We have found 
the same geometrical configurations discovered in [l[ . In 
the normal (nonphantom) case all solutions, except the 
de Sitter one, realize a three-dimensional spheroidal ge- 
ometry because the radial coordinate achieves a maxi- 
mum value rp, dubbed the equator. After the equator 



crossing, the same three scenarios studied in [l| may take 
place: a regular spheroid, a truncated spheroid having a 
scalar curvature singularity at a finite radial distance, 
and a closed spheroid having a Schwarzschild type singu- 
larity at the south pole, namely for \ = 7r - The presence 
of f3 does not give rise to any new scenario but slightly 
modifies the pressure and mass profiles in the three sce- 
narios mentioned above, as we have shown by studying 
their asymptotic behavior. 

We have also considered the possibility of a phantom 
gCg. In this case there is no equator and all solutions 
have the geometry of a truncated spheroid with the same 
type of singularity; namely the pressure diverges at a 
certain finite radius. We point out that the asymptotic 
behavior of the pressure near the resulting curvature sin- 
gularity is completely independent of (3. There are only 
two spherically symmetric static configurations in the 
presence of the phantom gCg: the first is a truncated 
spheroid, whereas the second is a wormholelike throat. 
The peculiar property of the latter is that it does not 
connect two asymptotically flat space-times. In these 
cases, the presence of (3 does not affect the qualitative 
behavior of the pressure. We have also shown that a 
smooth transition from the generalized to the standard 
Chaplygin gas near the asymptotic value p = —p = —A, 
required to preserve causality, does not affect the validity 
of all the previously obtained results. 

Finally, let us briefly recapitulate the main similari- 
ties and differences between the cases of the standard 
Chaplygin gas and of the gCg. First, the main result of 
our present investigation is that all geometrical configu- 
rations (and only these ones) which we have found in [l[ 
for the standard case are also present for the whole range 
of the gCg parameter < a < oo. Second, in spite of the 
qualitative similarity between the solutions for the gCg 
and the standard Chaplygin gas, there exists quantitative 
difference between these cases, especially for the normal 
(nonphantom) gCg, case 2 (see Sec. IIV|) . where the a de- 
pendence of the pressure profile and the geometric char- 
acteristics acquire a rather involved character. Third, 
another interesting aspect has been found in the case of 
the superluminal normal gCg. As we showed in for 
this case there is a good agreement with observation not 
only at the level of global cosmological dynamics, but also 
for the description of perturbations responsible for large 
scale structure formation. We have shown here that the 
technique employed in [9( to solve the problem of causal- 
ity in the superluminal fluid could be applied also in the 
context of the study of static spherically symmetric so- 
lutions of the TOV equations. 

A natural development of the techniques and results 
contained in this paper would be the analysis of a starlike 
object in the standard dynamical gCg universe and of a 
local spherical collapse of a gCg sphere, the latter process 
being a kind of generalization of the Lemaitre-Tolman- 
Bondi solution [H, H3 • Though we are fully aware of 
the difficulties of these problems, we are still hopeful to 
get some results in this direction in the near future. 
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